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The study of electron motion in semiconductor billiards has elucidated our understanding of 
quantum interference and quantum chaos. The central assumption is that ionized donors generate 
only minor perturbations to the electron trajectories, which are determined by scattering from 
billiard walls. We use magnetoconductance fluctuations as a probe of the quantum interference 
and show that these fluctuations change radically when the scattering landscape is modified by 
thermally-induced charge displacement between donor sites. Our results challenge the accepted 
understanding of quantum interference effects in nanostructures. 



PACS numbers: 72.20.My, 05.45.Mt, 73.63. Kv 

I. INTRODUCTION 

The use of solid-state devices to replicate elec- 
tron interference effects originally observed in vacuum 
holds considerable interest for fundamental and applied 
physics. The Aharonov-Bohm (AB) effect serves as an 
example [HE]- If a ring-shaped device enclosing an area 
A is connected on opposing sides to source and drain 
leads, then an electron flowing into the ring can take 
the clockwise or anticlockwise arm to exit at the other 
side. At temperatures low enough to establish electron 
phase coherence, quantum interference between these 
paths leads to a sinusoidal oscillation in conductance as a 
function of magnetic field B applied perpendicular to the 
ring with a period AB = h/eA. This geometry-induced 
AB effect was first observed using gold rings Sub- 
sequent measurements of semiconductor wires [IHZ] fo- 
cused on AB effects generated by material-induced scat- 
tering events. In these diffusive devices, a multitude of 
AB 'loops' arise from the many possible paths an elec- 
tron can follow through the device due to scattering both 
from the device walls and from electrostatic perturba- 
tions induced by ionized donors. The resulting magneto- 
conductance fluctuations (MCF) have a spectral content 
[7J determined by the distribution of loop areas and can 
therefore be used as a magnetofingerprint of the electron 
dynamics Jj. Semi-classical theory developed for diffu- 
sive devices emphasizes the crucial role played by donor 
perturbations by predicting that the magnetofingerprint 
can be changed significantly by changing the location of 
just one scattering site [SJO]- Experiments performed 
on diffusive wires have since confirmed the sensitivity of 
MCF to the movement of these scatterers [TdlfT5] . 

Ultrapure semiconductor heterostructures employ the 



modulation doping technique [16 to spatially sepa- 
rate the donors from the two-dimensional electron gas 
(2DEG) that forms at the heterostructure interface. This 
reduces the magnitude of the donor perturbations in- 
duced in the 2DEG's electrostatic landscape. Whereas 
the strong electrostatic perturbations induced by 'in- 
plane' donors in diffusive devices generate large-angle 
scattering events, these weaker perturbations are more 
likely to generate small-angle scattering. Because elec- 
tron mobility is naturally weighted towards large angle- 
scattering [IT] [18] , the associated mean free path can 
exceed the device size L, satisfying the traditional defini- 
tion of a semiconductor 'billiard' (L < i^) [T5]. However, 
this definition does not exclude the possibility of small- 
angle scattering within the device. It is not therefore a 
good gauge of ballistic transport, which requires electrons 
to follow straight trajectories between specular reflections 
from the device walls. This is highlighted by recent scan- 
ning gate microscopy (SGM) studies [I7J 1201 121] which 
show that small-angle scattering produces branching of 
trajectories at length scales much shorter than the mean 
free path for large-angle scattering. In particular, SGM 
performed on heterostructures with electron mobilities in 
excess of those used for fabricating billiards indicate that 
significant small-angle branching can be expected within 
semiconductor billiards |17j . 

These SGM studies re-ignite a central question for 
quantum transport research — how does the presence of 
small-angle scattering impact the AB effect in semicon- 
ductor billiards? The semi-classical theory developed for 
the diffusive regime does not quantify the minimum per- 
turbation size required for its applicability [9]. Conse- 
quently, it is not clear whether the sensitivity of the AB 
effect to perturbations highlighted by theoretical [5] and 
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experimental [TTHT5] investigations of diffusive devices 
extends to the smaller perturbations experienced by bil- 
liards. The pioneering experimental [22H26] and theo- 
retical [37] investigations of billiards assumed that donor 
perturbations did not significantly affect the electron dy- 
namics established by the scattering imposed by the de- 
vice walls. The experimental results presented here chal- 
lenge this picture by demonstrating that donor-induced 
small-angle scattering is the dominant factor influencing 
the MCF. 

We employ two techniques to relocate charge between 
donor sites — thermal activation and illumination with a 
LED. Thermal activation [TJ1 US] and illumination 
[l"3l l30l 131] have been used previously to demonstrate 
that charge relocation induces changes in the MCF for 
diffusive devices dominated by large-angle scatterers. In 
addition, illumination is a well-established technique for 
studying charge activation energies for donors [351 155] . 
Here, we apply these techniques to investigate the im- 
pact of charge relocation on the billiard's MCF. We show 
that redistribution of the ionized donors reconfigures the 
MCF despite the billiard's fabricated geometry remaining 
unchanged. This reappraisal highlights the challenges of 
realizing genuine quantum ballistic transport in semicon- 
ductor nanostructures. Furthermore, the results address 
fundamental aspects of quantum interference and its sen- 
sitivity to electron scattering [OJ . The MCF are observed 
to gradually reconfigure as the number of altered scat- 
tering sites in the billiard increases, demonstrating the 
cumulative impact of changing disorder on quantum in- 
terference processes. 

Our results also have broader implications beyond 
semiconductor physics. The quantum behavior of semi- 
conductor billiards can be used to investigate whether 
universal characteristics are evident across diverse sys- 
tems. For example, in quantum chaos studies, MCF were 
used to investigate the change from stable to chaotic elec- 
tron dynamics for billiards with differently shaped walls 
[2"2l |2"5] and for shape transitions within a single billiard 
[2H[5S]. These studies, which neglected donor scatter- 
ing, arc influential in the broader research community 
due to analogous phenomena in quantum wells [33]; in 
microwave [35] , optical [36j [37] , and acoustic [38] cavi- 
ties; and cold atoms in optical traps [39] , 
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FIG. 1. (Color online) Schematics of (a) GaAs/(AlGa)As 
and (b) (GaIn)As/InP billiards. The 2DEG (white) is lo- 
cated near the lower GaAs/(AlGa)As interface for (a) and 
in the (Gain) As layer for (b). The modulation-doping (red) 
consists of a 40 nm (AlGa)As layer with a Si donor density of 
1.2xl0 18 cm -3 for (a) and a 1 nm InP layer with a Si den- 
sity of 5xl0 18 cm" 3 for (b). The billiard walls (light blue) 
are denned by surface gates (yellow) in (a) and by etched 
trenches in (b) . The heterostructure in (b) has a uniform top 
gate separated from the surface by a 1 fim insulating layer 
(not shown), (c) Simulation of the energy band structure for 
the GaAs/(AlGa)As heterostructure, plotting the conduction 
band edge Ec v. distance from the 2DEG. (d) Magnetocon- 
ductance traces for a 1 /im (GaIn)As/InP billiard measured 
at temperatures T of (top to bottom) 12 K, 5 K, 1 K, and 
240 mK. The traces are offset for clarity. Inset: A plot of the 
electron phase coherence length 1$ as a function of T, and 
associated fit line, for the (GaIn)As/InP billiard. 



II. EXPERIMENTS 

To investigate the generic properties of electron in- 
terference in semiconductor billiards, we consider two 
widely studied material systems and two distinct meth- 
ods for defining billiards. In each CcLS6, cL 2DEG 
forms within a heterostructure (Fig. l(a, b)). For 
the GaAs/Al .33Ga .67As system, electrostatic deple- 
tion generated by patterned surface gates defines 'soft' 
billiard walls in the 2DEG (Fig. 1(a)) 00]. For 



the Gao.25Ino.75As/InP system, wet etching generates 
'harder' walls with steeper electrostatic gradients (Fig. 
1(b)) 01]. The 2DEG densities are n = 2 x 10 15 mT 2 
and 7 x 10 15 m~ 2 , giving Fermi energies of Ep = 8.4 
meV and 50 meV, respectively. Accordingly, the Fermi 
wavelengths (50 nm and 30 nm) are significantly smaller 
than the billiard widths of 1 /im, as required for semi- 
classical interference phenomena [27]. The fabrication 
of high-quality billiards involves the modulation doping 
technique shown in Fig. 1(c). Silicon donors are sepa- 
rated from the 2DEG plane by a spacer layer in order to 
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minimize the electrostatic perturbations induced in the 
2DEG by the ionized donors [16j 02] ■ Small-angle scat- 
tering induced by these perturbations is not sufficient to 
reduce to below the billiard width. The values for 
our billiards (3 /im and 6 /im for the GaAs/(AlGa)As and 
(GaIn)As/InP systems) are comparable to those mea- 
sured in the pioneering MCF experiments [22J [23] . 

Four-terminal conduction measurements were per- 
formed using AC lock-in techniques at a frequency of 37 
Hz and an excitation current of 1 nA. The bottom trace 
of Fig. 1(d) shows the MCF generated by wave inter- 
ference, superimposed on a classical conductance back- 
ground for a square (GaIn)As/InP billiard measured at 
T = 240 mK (base temperature) . Billiards support a dis- 
tribution of AB loop areas, leading to the observed range 
of frequencies in the magnetoconductance [221 126"! [2TI 141) . 
Raising T reduces the electron phase coherence length 
gU S3]: by T = 12 K, is smaller than the billiard 
width, suppressing the MCF and leaving only the classi- 
cal background. 

To investigate the impact of the donors on the AB 
loops, the billiards were cooled to base temperature in 
the dark and then warmed to temperatures high enough 
to supply the thermal energy necessary to relocate elec- 
tron charge between dopant sites. Within the 1 /im 2 
area of the (GaIn)As/InP and GaAs/(AlGa)As billiards, 
there are approximately 30,000 and 50,000 of these donor 
sites, respectively. Note, however, that these donors 
are distributed at various heights in the 40 nm modu- 
lation doped layer for the GaAs/(AlGa)As heterostruc- 
ture. This introduces considerable spatial overlap of the 
donors such that the associated electrostatic perturba- 
tions defined within the 2DEG generate a much smaller 
number of individual scattering events sites within the 
billiard. Simulations of GaAs/(AlGa)As heterostructures 
similar to ours indicate that the electrostatic landscape 
comprises collective features on the scale of 0.1 /.im, in- 
dicating approximately 100 scattering sites within the 
billiard [32]. If the trajectories forming the AB loops are 
determined purely by the fabricated geometry of the bil- 
liard walls, as traditionally believed, then the MCF will 
be immune to this shift in charge between dopants. 

We annealed the billiard devices by warming in the 
dark to an intermediate temperature Tj for a fixed time 
t = 30 minutes and then cooled the billiards back to base 
temperature. We then compared the pair of MCF traces 
taken at base temperature before and after this thermal 
cycle. Figure 2 shows three pairs of MCF traces mea- 
sured on the (GaIn)As/InP billiard. The bottom pair of 
traces were recorded back-to-back (i.e., Tj = 240 mK) 
and demonstrate the well-known reproducibility of MCF 
traces at low temperatures. The middle pair of traces 
correspond to Tj = 115 K and reveal the same high de- 
gree of correlation as the bottom pair. However, the 
top pair have visibly decorrelated after thermally cycling 
to Tj = 300 K. Our equivalent MCF measurements per- 



formed on the GaAs/(AlGa)As billiard, shown in Fig. 3, 
demonstrate that this decorrelation is a generic effect for 
semiconductor billiards. We note that the MCF of the 
GaAs/(AlGa)As billiard have a smaller high-frequency 
component than the (GaIn)As/InP billiard due to the 
billiard's smaller of 3.6 fim at T = 240 mK gj. 

To quantify the change in MCF as a function of Tj, we 
applied the following correlation function [TJJ [2SJ to each 
pair of traces G\(B) and G2{B): 
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[G 1 {B)-G 2 {B)\ 
N 



N = ([G X (B) - G y (B)} 



— , where (1) 
(2) 



The symbol ( ) b represents an average over the data 
points spanning — B c < B < B c , where B c is the field at 
which the cyclotron diameter matches the billiard width. 
For the GaAs/(AlGa)As billiard, B c = 0.16 T, and 640 
data points fall in this range; for the (GaIn)As/InP bil- 
liard, B c = 0.28 T, and 1120 data points fall in this range. 
The normalization constant N is calculated by averaging 
the correlations of 15 pairs of traces (G X (B) and G y (B)) 
that have been thermally cycled to Tj = 300 K. Adopting 
this normalization, the correlation scale varies between 1 
for mathematically identical traces to for decorrelation 
induced by a Tj = 300 K cycle. We note that the form of 
the correlation function minimizes noise in F for values 
close to 1. This is important for obtaining an accurate 
fit to the dependence of F on Tj (see below). 

In Fig. 4, we show the results for the (GaIn)As/InP 
and GaAs/(AlGa)As billiards. For comparison, we show 
previously reported decorrelations of the MCF for a 
GaAs wire [TJJ [25] ■ The wire was 10 /jm long, 50 nm 
high, and 90 nm wide, and was heavily doped with sil- 
icon (n = 5 x 10 24 m~ 3 , E F = 128 meV) distributed 
uniformly throughout the wire's cross-section. Due to 
the absence of modulation doping, the strong electro- 
static perturbations of these 'in-plane' dopants generates 
large-angle scattering, resulting in = 0.3 /im. The 
thermally-induced decorrelation of the MCF for this dif- 
fusive scattering system (strong perturbations) is clearly 
similar to that of the two modulation-doped billiards 
(weak perturbations). 



III. DISCUSSION 

For all three systems, we attribute the MCF decorre- 
lation to the relocation of electron charge among silicon 
dopants. This charge relocation changes the electrostatic 
landscape experienced by the electrons g2[ SJ] , which in 
turn reshapes the distribution of AB loops. The probabil- 
ity P for this process depends upon the available thermal 
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FIG. 2. (Color online) Comparison of MCF for the (Gain) As/InP billiard after being warmed to three intermediate temperatures 
Ti. Top to bottom, the Ti values for the pairs of traces are 300 K, 115 K, and 240 mK. The traces are offset for clarity. Top: 
Magnified comparison of traces taken with T = 115 K (lower [black] and middle [green]) and T = 300 K (lower [black] and 
upper [red]). 
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FIG. 3. (Color online) Comparison of MCF for the GaAs/(AlGa)As billiard after being warmed to three intermediate temper- 
atures Ti. Top to bottom, the T values for the pairs of traces are 300 K, 115 K, and 240 mK. The traces are offset for clarity. 
Top: Magnified comparison of traces taken with Ti = 115 K (lower [black] and middle [green]) and Ti = 300 K (lower [black] 
and upper [red]). 



energy according to P ~ exp(— /3/fcsTi), where j3 repre- billiard is held at Ti. We model the decorrelation be- 
sents the activation energy for charge relocation. The havior using F = exp[— ryexp(— /3/fc^Ti)], where rj — t/r a 
charge redistribution is also governed by the time t the and r is a characteristic time for charge transfer |28j . 
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The decorrelation data are fitted with this formula in 
Fig. 4. The insets in Fig. 4 are simulations demonstrat- 
ing how the form of the fall-off depends upon /? and -q: f3 
(upper) and rj (lower) are varied while keeping the other 
parameter fixed. The GaAs/(AlGa)As billiard's decor- 
relation commences at a higher Tj than for the other 
systems, indicating a higher j3 value. 

To explain this observation, consider the 
GaAs/(AlGa)As conduction band diagram of Fig. 
1(c) in more detail. Electrons undergo thermally- 
assisted tunneling from the 2DEG through the energy 
barrier at the GaAs/(AlGa)As interface to the doped 
layer's conduction band edge. In this process, a thermal 
energy of 130 meV is required for the electrons to reach 
barrier heights at which the barrier width is sufficiently 
narrow for tunneling. We note that the 20 nm tunneling 
distances required are similar to those observed in other 
tunneling structures |45j . Having tunneled into the 
(AlGa)As conduction band, DX centers [46] associated 
with the silicon dopants can then capture pairs of these 
electrons and change their charge state [J7]. This DX 
capture energy is ^250 meV for Alo.33Gao.67As [55] . 
This sequential process provides a possible explanation 
for the high /3 value of 350 ± 100 meV obtained from the 
fit. In contrast, in the (GaIn)As/InP heterostructure, 
the doped layer's conduction band edge resides only 
75 meV above the 2DEG [48 . Furthermore, the Si 
dopants are expected to form shallow traps 5.6 meV 
below the band edge [4*51. These smaller activation 
energies are consistent with the lower j3 value observed 
for this system (45 ± 20 meV). Finally, the (3 value of 
80 ± 20 meV observed for the GaAs wire is close to 
the expected activation energy of 100 meV, which is 
set by the height of the DX center above Ep for this 
layer [32]. We plan future MCF experiments to study 
the relocation of charge between donors and its detailed 
dependence on the entire thermal annealing cycle. The 
current experiments do not take into account the finite 
time to warm to X^. Future investigations will take into 
account activation of charge not just at Ti but also at 
the temperatures experienced during the warm-up and 
cool-down, which would allow more accurate 77 values to 
be measured. 

To confirm the dependence of MCF decorrelation on 
the transfer of electrons between the 2DEG and donor 
sites, we employed a LED as follows [32] . After cooling 
to base temperature in the dark, the (GaIn)As/InP bil- 
liard was illuminated using a series of LED pulses. The 
resulting transfer of electrons from the dopants to the 
2DEG was monitored with Hall measurements, which 
charted the rise in 2DEG density n from the initial value 
of n m i n to the saturation value of n max . The normal- 
ized n value, n norm = (n-n min )/(n max -n min ), was then 
measured in the dark as a function of rising temperature. 
The decrease in n norm shown in Fig. 5 results from the 
thermally-activated transfer of electrons from the 2DEG 
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FIG. 4. (Color online) A plot of F v. Tj for the GaAs wire (red 
circles), the (GaIn)As/InP billiard (light blue squares), and 
the GaAs/(AlGa)As billiard (green diamonds). The insets 
show the annealing function's form as f3 and r] are varied (see 
text). 
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FIG. 5. (Color online) A plot of F (left axis) v. T t (bot- 
tom axis) for the (GaIn)As/InP billiard (light blue squares), 
as well as n llor m (right axis) v. T (top axis) for the 
(GaIn)As/InP billiard (blue crosses). The decorrelation data 
for the (GaIn)As/InP are identical to those shown in Fig. 4. 



back to the donors. The form of n norm v. T is strikingly 
similar to that of F v. Tj, highlighting the dependence of 
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MCF decorrelation on charge transfer from the 2DEG. 

We note that whether electrons have to tunnel from 
the 2DEG to the (AlGa)As layer in order to be captured 
by the donor sites, as observed here, depends on the pre- 
cise layer architecture of the heterostructure. For ex- 
ample, for heterostructures with much wider (AlGa)As 
layers than those employed for billiards research, par- 
allel conduction can occur in the (AlGa)As layer [50] . 
Consequently, for these heterostructures, electron trans- 
fer across the interface barrier is not necessary to popu- 
late the (AlGa)As layer and the donors [50"] . 

IV. CONCLUSIONS 

In conclusion, we have shown that the electrostatic 
profile of ionized donors strongly affects the interference 
of electron waves in semiconductor billiards. By demon- 
strating this effect for two material systems and different 
fabrication methods, we propose this is a generic limita- 
tion of semiconductor billiards. This contradicts previ- 
ous models in which the quantum dynamics was inter- 
preted as being determined purely by the billiard walls 
[22H27] EJJ and which predict that the MCF should be 
robust to thermal cycling. In contrast to this traditional 
billiard picture, our experiments show that the modu- 
lation doping technique does not reduce the donor elec- 
trostatic perturbations sufficiently to ignore their pres- 
ence. In particular, we have demonstrated experimen- 
tally that the sensitivity of the AB effect, highlighted by 
semi-classical theory developed for diffusive regime [3], 
extends to billiards even though the size of the electro- 
static perturbations in billiards is reduced significantly 
compared to diffusive systems. 

In addition to revealing the role of donor-induced dis- 
order, our results also highlight central aspects of the 
sensitivity of quantum interference to electron scatter- 
ing. The semi-classical theory predicts that changing 
one scattering site can have a significant effect on the 
MCF [S] , raising the question of whether one site change 
in the billiard is sufficient to completely decorrelate the 
MCF. The observed gradual decorrelation indicates that 
the total number of altered sites is important and has a 
cumulative effect on the MCF. 

Our experiments can be combined with earlier obser- 
vations demonstrating the lack of MCF sensitivity to bil- 
liard wall properties such as geometry and the softness of 
the wall's electrostatic profile [41] [52]. Taken together, 
these results can be seen as supporting an earlier pro- 
posal that the electron dynamics in billiards are gener- 
ated by donor scattering [IT]. Within this picture, the 
walls serve to repeatedly reflect electron trajectories to- 
wards the donor scattering sites [UJ. However, we note 
that some MCF features, such as those induced by sym- 
metries in the billiard walls 51j, nevertheless display ge- 
ometry dependences. Our results therefore highlight the 



need for quantitative investigations aimed at determining 
the relative contributions of walls and ionized donors to 
the MCF, which might be influenced by trajectory length 
and stability [S3], as well as dynamical effects such as 
quantum scarring [53] . 

In contrast to the sensitivity of the MCF to donor dis- 
tributions, the classical magneto conductance is expected 
to be determined predominantly by the billiard geome- 
try. This difference in conduction processes originates 
from the sensitivity of the AB effect. The magnetic field 
dependence of the classical conduction is due to 'focus- 
ing' of the trajectories into the billiard's entrance and exit 
due to curvature from the Lorentz force [19]. Classical 
focusing is dependent only upon the opening by which 
the trajectories leave the billiard and is therefore rela- 
tively insensitive to small-angle scattering events. The 
AB effect is instead sensitive to the precise routes taken 
by the electrons through the billiard. Thus, whereas the 
presence of donors in the billiard will simply perturb the 
classical condition by broadening the focusing features, 
the donors will change the scattering loops substantially. 
This insensitivity to donors explains the previous success 
in observing classical ballistic effects [TJ0 [55] [5H] . 

Our observations demonstrating the importance of 
donors for quantum transport through billiards is dif- 
ferent from earlier results performed on quantum point 
contacts, where the sensitivity to donors can be ex- 
plained in terms of one-dimensional channels resulting 
from quantum confinement [571460] . These results on 
QPCs, along with our results on billiards, emphasize the 
central role that disorder plays for quantum transport 
in solid-state environments: the associated lack of repro- 
ducibility between devices and thermal cycles has impli- 
cations not only for fundamental research but also for 
future nanoscale electronic devices. We expect this to 
apply to novel systems such as graphene 161,1621 and car- 
bon nanotubes [63]. 
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